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Abstract 

The strong electromagnetic fields in peripheral heavy ion collisions 
give rise to photon-photon and photon-nucleus interactions. I present 
a general survey of the photon-photon and photon-hadron physics ac- 
cessible in these collisions. Among these processes I discuss the nuclear 
fragmentation through the excitation of giant resonances, the Coulomb 
dissociation method for application in nuclear astrophysics, and the 
production of particles. 



1 Peripheral Heavy Ion Collisions 

The held of peripheral atomic collisions was born in 1924, when E. Fermi 
had the ingenious idea of relating the atomic processes induced by fast 
charged particles to processes induced by electromagnetic waves. In 1934- 
1935, Weizsacker and Williams modified Fermi's calculation by including the 
appropriate relativistic corrections. The original Fermi's idea is now known 
as the Weizsacker- Williams method [2, 3], an approximation widely used 
in coherent processes in atomic, nuclear and particle physics. In Fermi's 
method the electromagnetic (strong) field generated by a fast particle is 
replaced by a flux of equivalent photons (flux mesons or gluons). 

In peripheral heavy ion collisions (PHIC) the number of equivalent pho- 
tons, n(u) of energy u can be calculated classically, or quantum-mechanically. 
For the electric dipole (El) multipolarity both results are identical under the 
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assumption of very forward scattering [4]. In ref. [4] the number of equiv- 
alent photons for all multipolarities was calculated exactly. It was shown 
that for the electric dipole multipolarity, El, the equivalent photon number, 
iiei(w), coincides with the one deduced by Weizsacker and Williams. It was 
also shown that in the extreme relativistic collisions the equivalent photon 
numbers for all multipolarities agree, i.e, uei{uj) ~ riE2(w) ~ «mi(^) ~ •••• 
This is shown in figure 1. 

According to Fermi's idea, the cross sections for one- and two-photon 
processes depicted in figure 2(a,b) are given by 

f n(co) 7 f n(ui)n(u 2 ) 77 , , 
&x = / duj a x {ui) , and ax = / «<^i«<^2 <?x (wi,^) , 

(1) 

where a x (u) is the photon-induced cross section for the energy u, and 
(Tx (^1,^2) is the two-photon cross section. Note that we do not refer to 
the photon momenta. The virtual photons are real: q 2 = 0, a relation always 
valid for PHIC. 

At relativistic bombarding energies one has (in natural units) 

where 7 = (1 — the relativistic Lorentz factor in the frame of 

reference of the target, Z p e is the projectile charge and a = 1/137. 

Eq. 1 is called the equivalent photon approximation (EPA) , or Weizsacker- 
Williams method. It is valid for processes c, d, e, and / in figure 2. But the 
approximation is not applicable to processes a and b in figure 2, in which 
multiple virtual photons are exchanged between the projectile, the target, 
and/or with the produced particle(s). 

For one-photon processes, e.g., Coulomb fragmentation, a x (u) is local- 
ized in a small energy interval and one gets a cross section in the form 
a = A In 7 + B, where A and B are coefficients depending on the system. 
The In 7 factor is due to the equivalent photon number, n(uj), which is 
approximately independent of to in the integral range of interest. For two- 
photon processes, besides the In 2 j c from n\ and 722, a third hi7 c arises from 
the integral over uj\ {lo\ and L02 are related by energy conservation). Note 
that here we used 7 C of a heavy ion collider, so that 7 = 27^ — 1, with 7 C 
the collider Lorentz gamma factor (e.g., 7 C ~ 100 for the RHIC collider at 
Brookhaven) . 
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Figure 1: Virtual photon numbers (in units of Z 2 a) for the El, E2 and Ml 
multipolarities. Upper (lower) figure is for 7 = 2 (7 = 100). to is the photon 
energy and R is the sum of the nuclear radii. 



Most applications of PHIC were reviewed in ref . [5] . Since then a great 
amount of work has been performed in this field. The coherent 7-7 and 7- A 
interactions in very peripheral collisions at relativistic ion colliders has been 
reviewed recently in ref. [6]. In this review article I will present a general 
account of the theory and the latest developments in the field of PHIC with 
fixed targets at intermediate energies ~ 100 MeV/nucleon), and at 

ultra-relativistic energies of present day (the Relativistic Heavy Ion Collider 
(RHIC), at Brookhaven) and future heavy ion collider (the Large Hadron 
Collider (LHC), at CERN) energies . 
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(a) „ (b) (c) 




Figure 2: PHIC processes: (a) one-photon, (b) two-photon, (c) 
Bremsstrahlung, (d) Delbriick scattering, (e) pair-production, and (f) pair- 
production with capture. 

2 Relativistic Coulomb Excitation and Fragmen- 
tation 

Relativistic Coulomb excitation is now a popular tool for the investigation of 
the intrinsic nuclear dynamics and structure of the colliding nuclei, specially- 
important in reactions involving radioactive nuclear beams [1, 7, 8, 9, 10, 11, 
12]. The advantage is that the Coulomb interaction is very well known. The 
disadvantage is that the contribution of the nuclear-induced processes also 
play a role in some situations. The treatment of the dissociation problem by 
nuclear forces is very model dependent, based on eikonal or multiple Glauber 
scattering approaches [7, 13, 14, 15, 16]. Among the uncertainties are the 
in-medium nucleon-nucleon cross sections, the multiple scattering process 
and the separation of stripping from elastic dissociation of the nuclei [16]. 
Nonetheless, specially for the very weakly-bound nuclei, relativistic Coulomb 
excitation has lead to very exciting new results [7, 8, 13, 14, 15]. 

2.1 Reactions with Radioactive Beams 

Coulomb breakup of weakly-bound nuclei may involve single or multiple 
photon-exchange between the projectile and the target. In the first case, 
perturbation theory gives a direct relation between the data and the matrix 
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Figure 3: S-factors (Sn) for the 7 Be(p, "f) 8 B reaction. The inset is a blowup 
of the low energy region. See text for details. 

elements of electromagnetic transitions. Such matrix elements are the clear- 
est probes of the nuclear structure of these nuclei, since one cannot perform 
experiments with real photons or with electron scattering off nuclei far from 
the stability valley. In the second case, often called by post-acceleration, 
or reacceleration, effects [13, 14, 15, 17, 18], one has to perform a non- 
perturbative treatment of the reaction what complicates the extraction of 
the electromagnetic (mainly El) matrix elements. 

Among several possible studies of interest, one expects to learn if the 
Coulomb-induced breakup proceeds via a resonance or by the direct dissoci- 
ation into continuum states [13, 14, 15]. There is a strong ongoing effort to 
use the relativistic Coulomb excitation technique also for studying bound ex- 
cited states in exotic nuclei, to obtain information on gamma-decay widths, 
angular momentum, parity, and other properties of hitherto unknown states 
[8, 11, 12]. 

Another application of peripheral collisions with radioactive nuclear beams 
is in astrophysics. Radiative capture reactions are known to play a major 
role in astrophysical sites, e.g., in a pre-supernova [19, 20]. Some reactions 
of interest for astrophysics, e.g., 7 Be (p, 7) 8 B, can be studied via the inverse 
photo-dissociation reaction 8 B(j,p) 7 Be [21]. Since the equivalent photon 
numbers in eq. 1 can be calculated theoretically , an experimental mea- 
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surement of the Coulomb breakup reaction a + A — ► b + c + A is useful 
to obtain the corresponding 7 induced cross section 7 + a — ► b + c. Using 
detailed balance, this cross section can be related to the radiative capture 
cross section b + c — ► a + 7, of astrophysical interest [21] (see eq. 4). 

The radiative capture cross sections in nuclear astrophysics are often 
written in terms of the astrophysical S-factor [19, 20], defined by 

S(E) = Eg (E) exp [2ttt] (E)] , (3) 

where r\ (E) = ZbZ c e 2 /h^/2[it }C E , and E is the relative kinetic energy of the 
nuclei b and c in the reaction b(c, 7)0. The cross section for the dissociation 
process a + A — ► b + c + A can be written as 



da ^ bc ° 2 <t? \ (2Jb + l)(2J c + l) , . r , . . 

IE = ~Ef n{Ex) (2XTl) 6XP h2?rr/ (E)] ' (4) 

where ^ x = E + B is the sum of the relative energy E and the binding 
energy B of the two fragments. Jj is the ground-state angular momentum 
of the nucleus i. 

The Coulomb dissociation method is more useful when higher order ef- 
fects are under control, so that the eq. 4, obtained in lst-order perturbation 
theory, is valid. Higher order effects can be taken into account in a coupled 
channels approach, or by using higher order perturbation theory. 

Expanding the nuclear wave function in the set {[ j); j = 1,N} of 
eigenstates of the intrinsic Hamiltonian Hq, where N is the number states 
included in the calculation, one obtains the coupled equations for the occu- 
pation amplitude of the state k as 



ih a k {t) = ( k I V(t) I j) exp [i(E k - Ej)t/h] aj(t) , k = 1 to N , 



N 

It ' 

(5) 

where E n is the energy of the state \n) . This set of equations can be solved 
numerically if one has a theoretical model or experimental information on 
the matrix elements (k \ V(t) \ j). In the most relevant situations the 
perturbing potential is the Coulomb interaction of the fragments with the 
target. This interaction is expanded into electric and magnetic multipolar- 
ities. One thus needs information on the electromagnetic matrix elements 
between the bound states and the continuum states of the nucleus to carry 
out the coupled-channels calculation. This approach has been used in refs. 
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Table 1: Radiative capture reactions of interest in nuclear astrophysics which 
can be studied with the Coulomb dissociation method. 



Reaction 


Ti/2 (projectile) 


Astrophysical application 


■ i He(a,j) Be 


53.3 days 


Solar-neutrinos 


7 Be (p, 7 ) s B 


770 ms 


3 TT 1 1 

°He abundance 


7 Be (a,j) LL C 


20.4 min 




A tt / T \b r • 

4 iJe (d, 7) Li 


Stable 


Primordial nucleosynthesis 


b Lz (p, 7) Be 


53.3 days 




fir-/ \ 1U t-» 

b Li(a,-f) B 


Stable 




4 He (t,7)' Li 


Stable 




7 t ■ / \ 11 7-1 

7 Li(a^) B 


Stable 




n B(p,7) C 


Stable 




9 Be (p, 7) B 


Stable 




in n / \11 n 


20.4 min 




7 Li (n, 7)* 


842 ms 


Inhomogcncous Big Bang 


8 Li(n,-y) Li 


178 ms 




12 C(n, 7 ) i3 C7 


Stable 




14 C(n, 7 ) is C7 


2.45 s 




14 C(a, 7 ) is O 


Stable 




12 C(p, 7 ) i3 ^V 


10 min 


CNO cycles 


16 0(p, 7 ) iY F 


65 s 




13 7V(p, 7 ) i4 


70.6 s 




20 7Ve(p,7)^iVa 


22.5 s 




U C(P,7)"JV 


11 ms 


Hot p-p chain 


15 0(a, 7 ) iy iVe 


17.2 s 


rp process 


31 S(p,7) 32 


291 ms 




12 C(a, 7 ) ib O 


Stable 


Helium burning 


16 0(a, 7 ) zu iVe 


Stable 




14 Ar(a, 7 ) is F 


109.7 min 




22 M 5 (p,7)^A/ 


3.86 s 


rp bottlenecks 
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[22, 23] to treat the problem of Coulomb reacceleration of fragments follow- 
ing the breakup. More recently, a similar technique has been used in ref. 
[24]. 

Another approach to the Coulomb postacceleration, or reacceleration, 
problem is to integrate the time-dependent Schrodinger equation directly 
for a given model Hamiltonian. This approach is only useful when one can 
use a two-body potential model for the nucleus. Then, expanding the two- 
body wavefunction into angular components one gets the time-dependent 
wave equation 



dr 2 



1(1 + 1) _ 2^c 

r 2 /j 



V{r) 



I'm' 



I x 2^ fec dui m 
l 'U Vm \r) = — (6) 



where u\ m (r) is the radial part of the Zm-component of the two-body wave- 
function. Sjj^, is a source term for the postacceleration, arising from the 
multipole component Ira of the Coulomb field of the target. This approach 
has been used in ref. [15] and later developed in refs. [17, 18, 25, 26, 27]. 

As an example of application of the Coulomb dissociation method, we 
show in figure 3 the result of an experiment performed at the GSI laboratory, 
in Darmstadt, Germany [28]. The S-factor obtained in this experiment is 
shown in figure 3 as solid circles, by a direct application of eq. 4, using 
Sn(E) obtained from theoretical models. The solid curve is a fit using a 
theoretical model of ref. [23] , whereas the dashed curve is a calculation done 
in ref. [29]. 

Alternative experiments have obtained S±j by measuring the parallel 
momentum distribution of the fragments after Coulomb breakup [33, 32]. 
This technique was applied to several experiments at the NSCL laboratory, 
of Michigan State University [32, 31, 30]. The technique also allows the 
disentanglement of the El and E2 contributions to the breakup by looking 
at the asymmetries in the momentum distribution following the Coulomb 
breakup of the projectile [33]. 

In table I a set of reactions are shown which can be studied with the 
Coulomb dissociation method [34]. 

Besides the reaction 7 Be (p, 7) 8 B, the Coulomb dissociation method has 
also been applied to the study of 4 He (d, 7) 6 Li (see ref. [35]), 12 C (n, 7) 13 C 



(see ref. [36, 37]), n C (p, 7) 12 N (see ref. [38]), and 12 C(a, 7 ) iD (see ref. 
[39]). 



,16 
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2.2 Multiphonon Resonances 

Giant dipole resonances (GDR) occur in nuclei at energies around 10-20 
MeV. Assuming that they are harmonic vibrations of protons against neu- 
trons, one expects that DGDRs (Double Giant Dipole Resonances), i.e., two 
giant dipole vibrations superimposed in one nucleus, will have exactly twice 
the energy of the GDR [5, 9, 10]. 

Assuming that one knows a^ DR (E) somehow (either from experiments, 
or from theory), a simple harmonic model can be formulated to obtain the 
Coulomb excitation cross sections of these states. In this model the inclusion 
of the coupling between all multiphonon states can be performed analyti- 
cally [40]. Following the same reasoning leading to eq. 1, to first-order the 
probability to excite a giant resonance state at energy E is given by 



p^ {Ejb) = ^R a GDR {Eh (?) 
E 



where N(E, b) is the number of equivalent photons per unit area in a collision 
at impact parameter b. The relation between n(E) of eq. 1 and N(E, b) is 
given by n(E) = 2vr / N(E, b) b db. 

In the harmonic oscillator model the excitation probabilities calculated 
to first-order, are modified to include the flux of probability to the other 
states [5]. That is, 

P(E, b) = P lst (E,b) exp{-P ls '(6)} , (8) 

where P lst (b) is the integral of P(E,b) over the excitation energy E. In 
general, the probability to reach a multiphonon state with the energy 
from the ground state, with energy E^°\ is obtained by an integral over all 
intermediate energies [41] 

p(») ft) = 1 exp{-P ls *(6)} J d£ (n_1) dE (n ~ 2) ... dE {1) (9) 

x plst( E (n) _ E {n-1)^ plst( E (n-l) _ E (n-2) ^ _ pistil) _ E (0) jfify 

Integrating this equation over impact parameters yields the cross section for 
the excitation of the multiphonon state n, as a function of the excitation 
energy E^ n \ 

A series of experiments at the GSI laboratory have obtained the energy 
spectra, cross sections, and angular distribution of fragments following the 



9 



decay of the DGDR [42, 43, 44, 45, 46, 47, 48]. It was shown that the 
experimental cross sections are about 30% bigger than the theoretical ones. 
This is shown in figure 4 where the cross sections for the excitation of 1- 
phonon (GDR), 



01 ~ 27rS'ln 
while for the 2-phonon state it is 



2 7 4 /3 
A 1/3 + A 1/3 



(11) 



o2 

a~0.1 5 , (12) 

(4 /3 +4 /3 ) 2 mb 

where 

g = 5.45xl0- 4 Z ^f r mb. (13) 

The dashed lines of figure 4 are the result of more elaborate calculations 
[9, 10, 49, 50, 51, 52, 53]. 

The GSI experiments are very promising for the studies of the nuclear 
response in very collective states. One should notice that after many years 
of study of the GDRs and other collective modes, the width of these states 
are still poorly explained theoretically, even with the best microscopic ap- 
proaches known sofar. The extension of these approaches to the study of 
the width of the DGDRs will be helpful to improve such models [10]. 

In heavy ion colliders the mutual Coulomb excitation of the ions (leading 
to their simultaneous fragmentation) is a useful tool for beam monitoring 
[54]. A recent measurement at RHIC [55], using the Zero Degree Calorime- 
ter to measure the neutron decay of the reaction products, has proved the 
feasibility of the method. The theoretical prediction of about 3 b for this 
process, agrees quite well with the experimental results. 

The DGDR contributes to only about 10% of the total fragmentation 
cross section induced by Coulomb excitation in PHIC. The main contribu- 
tion arises from the excitation of a single GDR, which decays mostly by 
neutron emission. This is also a potential source of beam loss in relativistic 
heavy ion colliders [1] , and an important fragmentation mode of relativistic 
nuclei in cosmic rays. 
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Figure 4: Cross sections for the excitation of the GDR and the DGDR for several 
systems. The dashed curves are theoretical calculations. See text for more details. 



3 Atomic Processes 



3.1 Atomic ionization 



The cross sections for atomic ionization in PHIC are very large, of order 
of kilobarns, increasing slowly with the logarithm of the RHI energy. For 
a fixed target experiment using naked projectiles one gets [5] (here we use 
k = c = l) 



aj = C(3) 



/j P r e 



1.8vr + 9.81n ( 

\Ztol ) 



(14) 



which decreases with the target charge Zt- This is due to the increase of 
the binding energy of if-electrons with the atomic charge. The probability 
to eject a .fT-electron is much larger than for other atomic orbitals. In 



this equation 



' /m is the classical electron radius, and £(3) = 1.202 



accounts for the ionization of electrons from higher orbits. 

The first term inside brackets in eq. 14 is due to close collisions assuming 
elastic scattering of the electron off the projectile, while the second part is 
for distant collisions, with impact parameter larger than the Bohr radius [5]. 
Recently, Baltz [56] has shown that the numerical factors in the equation 
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Figure 5: Atomic ionization cross sections for Pb 81+ (33 TeV) projectiles 
impinging on several targets. The solid and dashed curves are theoretical 
calculations. See text for more details. 

above should be replaced by 1.8 — > 1.74 — 1.83 and 9.8 — > 7.21, respectively, 
when one treats the electronic wavefunctions with the relativistic corrections. 

Recent experiments have reported ionization cross sections for Pb 81+ 
(33 TeV) beams on several targets [57]. In this case, the role of projectile 
and target are exchanged in the previous equation. In figure 5 show the 
results of eq. 14 (dashed line) are compared to the experimental data. Since 
the targets are screened by their electrons, the discrepancy is expected. 
Calculations by Anholt, Becker and collaborators [58, 59, 60] (solid line) or 
of Baltz [56] also yield larger cross sections than the experimental data. 

Non-perturbative calculations, solving the time dependent Dirac equa- 
tion exactly, were first performed by Giessen and Oak Ridge groups [61, 62]. 
The main problem is to adequately treat the several channels competing 
with the ionization process, specially for atoms with more than one electron. 
Also, the effects of screening (static and dynamic) are hard to calculate. On 
the experimental side, there are little data available for a meaningful com- 
parison with theory. However, atomic ionization should be taken seriously 
in accelerator design when a residual gas remains in the accelerator pipes. 
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3.2 Bremsstrahlung and Delbriick scattering 



Bremsstrahlung (fig. 2c) is a minor effect in PHIC. It has been calculated 
in ref. [5, 63]. The cross section is proportional to the inverse of the square 
mass of the ions. Most virtual photons have very low energies. For 10 MeV 
photons the central collisions deliver 10 6 more photons than the peripheral 
ones [63]. 

For a collider the Bremsstrahlung differential cross section is given by 



where m^v is the nucleon mass, 7 = 27^ — 1, where 7 C is the collider Lorentz 
gamma factor (j c ~ 100 for RHIC/BNL), and R is the nuclear dimension 
(R ~ 2.4 x A 1 / 3 fm) [63]. 

For very low energy photons (oj ~ 100 eV) the whole set of particles in 
a beam bunch can act coherently and a great number of Bremsstrahlung 
photons can be produced. This has been proposed as a tool for monitoring 
the bunch dimension in colliders [64]. Recently, an experiment has been 
approved at RHIC/BNL to measure this kind of coherent effect, including 
the emission of real photons by the interaction of a bunch with the edges of 
the accelerator magnets [65]. 

Delbruck scattering (7* + 7* — > 7 + 7) involves an additional a 2 as 
compared to pair production and has never been possible to study experi- 
mentally. The cross section is about 50 b for the LHC [63] and the process 
is dominated by high-energy photons, E 1 S> m e . A study of this process 
in PHIC is thus promising if the background problems arising from central 
collisions can be eliminated. No experiments of Bremsstrahlung or Delbruck 
scattering in PHIC have been performed so far. The total cross section for 
Delbruck scattering (to 3> m e ) in colliders is given by [63] 



4 The Electromagnetic Field of an U It rarelat i vis- 
tic Charge 

As discussed above, the calculation of processes involving the multiple ex- 
change of photons is not as easy to perform. For example, substantial dif- 
ferences have been found between the cross sections for the production of 



d<7 7 
dto 




) 



(15) 




(16) 
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e + e~ pairs calculated within several approaches [67, 68]. However, the cal- 
culation of multiple photon exchange can be considerably simplified at very 
high bombarding energies if one assumes that the electromagnetic field gen- 
erated by the projectile (or target, depending on the reference frame) lies 
entirely on a plane perpendicular to its trajectory [69]. We show next that 
this approximation is well justified [70]. 

The scalar part of the electromagnetic potential generated by an ultra- 
relativistic particle is proportional to V (p, z,i) = (f>( p, z, t) where 4> (p, z, t) 
is the Lienard-Wiechert potential at a point r = (p, z), generated by a rela- 
tivists particle with velocity v = vz and impact parameter b (in our units 
H = c = m e = 1), 

(p, z, t) = jZa [(b - p) 2 + 7 2 (z - vtf] ~ 1/2 . (17) 

In these expressions 7 = (l — v 2 ) 1//2 is the Lorentz contraction factor. The 
vector part of the electromagnetic interaction is proportional to eq. 17 and 
, for simplicity, we leave it out of this proof. 

Using the Bethe-integral [5], the potential (17) can be written in the 
form 

1 f g-iq-Ugiq.r 

4> (p, z,t)=Za^J d\ q2 _ y2ql , (18) 

where u = b + vi and q = (q<, q z ). The integral in (18) diverges logarith- 
mically as v — > 1. 

For practical purposes one needs the interaction in the form V (p, z, t) = 
4>(p,z,t) -cf)(0,z,t), i.e., 

V (p, z,t)=Za— 2 j d\ e^ L^ • (19) 
Let us define 

V(p,z,t) 



Za 
If 1 

= - d 2 q t — exp(-zq t .b) [exp(iq t .p) - 1] $ z (q t , z, t)(20) 
vr J qf 



Then 



2 r iq z (z-vt) 

*.(*. m) = I f^^^.m —1— •. (2D 
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Now, using lirriA^+oo (A/2) e~ A ^ = 5 (x), we see that for 7 — > 00, <E> 2 does not 
depend on q t , and assumes the form of a delta function: & z (z, t) = 5 (z — vt). 
In this limit, we can write (20) as 



*(p,z,t) =6(z-t)* p {p,t) , (22) 



with 



If 1 

$p (P, t) = - d 2 q t -3 exp (-zq*.b) [exp (iq t .p) - 1] 
^ J Qt 

= 2 / ^{J [ gt | p -6|]_ J ( % 6)} , (23) 



where Jo is the cylindrical Bessel function. The integral over each Bessel 
function diverges, but their difference does not. To show this one regularizes 
the integrals by using 



, xJn (ax) T , . , . 

dX 2 , 2 = K ( ak ) . ( 24 ) 



+ fc 2 



where i^o is the modified cylindrical Bessel function. Taking the limit k — > 0, 
and using (afc) ~ In (ak), for small values of ak, one gets 

$(p,t)=ln^^£. (25) 

This is the solution of the Coulomb potential of a unit charge in 2- 
dimensions. An easy way to see this is to use Gauss law for the electric field 
in two dimensions. One obtains E ~ 1/6, where 6 is the distance to the 
charge. Since E = —dQ/db, the logarithmic form of $ is evident. 

The above derivation illustrates the validity of the approximation in 
terms of the transverse momentum transfer qt. It should fail for very soft 
processes, i.e., those for which q t — > 0. Also, it requires that q z is small 
compared to 7. As shown in ref. [5], q z values up to 7 units of the electron 
mass contribute appreciably to the integrals involved in the production of 
free, and of bound-free, e + e~ pairs. 

The delta-function interaction yields reasonable results as long as ujb/j < 
0.1. [70]. This amounts to 6 < O.Ty/a;. As shown in ref. [5], the most 
effective impact parameters for this process are of the order of 6 ~ 1/m. In 
ref. [70] it has been shown that the differential cross sections for e + e~ pair 
production are well described up to energies of the order of 0.17m, by using 
the delta-function interaction. 
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For other situations, e.g., nuclear fragmentation due to the electromag- 
netic interaction, the most effective impact parameter is given by b ~ R, 
where R ~ 10 fm. We thus expect that the delta-function interaction works 
well for e < O.I7 MeV. Note that 7 is huge for RHIC and LHC energies, and 
the approximation works well for all energies of practical interest in nuclear 
fragmentation. 

The basic idea of the delta-function interaction is that the electromag- 
netic field of a relativistic charge looks like a very thin pancake. Those 
processes which do not involve too large energy transfers, will not be sensi- 
tive to the spatial variation of the field. Then the delta-function is a good 
approximation. For typical energy transfers of the order of 10-100 MeV in 
nuclear fragmentation, the approximation works well for b < OAj/u fm. To 
calculate total cross sections it is always necessary to account for those large 
impact parameters at which the delta- function approximation fails. Similar 
conclusions have been drawn in a recent article on projectile-electron loss in 
relativistic collisions with atomic targets [71]. 

4.1 Free and bound-free electron-positron pair production 

It has been demonstrated long time ago [72, 73, 74, 75] that to leading order 
in In 7, the e + e~-pair production in PHIC is given by 

-=^^e 2 ln 3 g) • (26) 

A renewed interest in this process appeared with the construction of 
relativistic heavy ion accelerators. For heavy ions with very large charge 
(e.g, lead, or uranium) the pair production probabilities and cross sections 
are very large. They cannot be treated to first order in perturbation theory 
[5], and are also difficult to calculate. This resulted in a great number of 
theoretical studies [76, 77, 78, 79, 80, 81]. 

Replacing the Lorentz compressed electromagnetic fields by delta func- 
tions, and working with light cone variables, one has developed more elabo- 
rate calculations [82], recently. One still debates upon the proper treatment 
of Coulomb distortion of the lepton wavefunctions, and of production of 
n-pairs [82]. 

An important phenomenon occurs when the electron is captured in an 
atomic orbit of the projectile, or of the target. In a collider this leads to 
beam losses each time a charge modified nucleus passes by a magnet down- 
stream [63] . A striking application of this process was the recent production 
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Figure 6: Pair production with capture for Pb S2+ (33 TeV) beams on several 
targets. The solid and dashed curves are theoretical calculations. See text 
for more details. 

of antihydrogen atoms using relativistic antiproton beams [83, 84]. Here 
the positron is produced and captured in an orbit of the antiproton. Early 
calculations for this process used perturbation theory [5, 58, 85, 86, 87, 88]. 
Some authors use non-perturbative approaches, e.g., coupled-channels calcu- 
lations [76, 77, 78, 79, 80, 81]. Initially some discrepancy with perturbative 
calculations were found, but later it was shown that non-perturbative cal- 
culations agree with the perturbative ones at the 1% level (see, e.g., ref. 
[82, 89, 90, 91, 92, 68, 93]). 
The expression 

° = exp(2vrZ«)-l T {^) ~ s\ (27) 

for pair production with electron capture in PHIC was obtained in ref. [5]. 
The term in the denominator is the main effect of the distortion of the 
positron wavefunction. It arises through the normalization of the contin- 
uum wavefunctions which accounts for the reduction of the magnitude of 
the positron wavefunction near the nucleus where the electron is localized 
(bound). Thus, the greater the Z, the less these wavefunctions overlap. Eq. 
27 predicts a dependence of the cross section in the form a = A In 7 + B, 



17 



where A and B are coefficients depending on the system. This dependence 
was used in the analysis of the experiment in ref. [57]. 

In recent calculations, attention was given to the correct treatment of 
the distortion effects in the positron wavefunction [95]. In figure 6 we show 
the recent experimental data of ref. [57] compared to eq. 27 and recent 
calculations [94, 95]. These calculations also predict a In 7 dependence but 
give larger cross sections than in ref. [5]. The comparison with the exper- 
imental data is not fair since atomic screening was not taken into account. 
When screening is present the cross sections will always be smaller up to a 
factor of 2 (see ref. [5]). The conclusion here is that pair production with 
electron capture is a process which is well treated in first order perturbation 
theory. However, eq. 27 is shown to be only a rough estimate of the cross 
section [94]. As with the production of free pairs, the main concern here is 
the correct treatment of distortion effects (multiphoton scattering) [94, 95]. 

The production of para-positronium in heavy ion colliders was calculated 
[66]. The cross section at RHIC is about 18 mb. This process is of interest 
due to the unusual large transparency of the parapositronium in thin metal 
layers. 

The production of mesons in two- and three-photon collisions can be 
studied following the same formalism as used in QED for the production of 
positronium [96]. This is shown in the next section. 

5 Photon-photon Fusion Mechanism 
5.1 Meson Production 

The production of heavy lepton pairs or r + r _ ), or of meson pairs 

(e.g., 7r + 7r~) can be calculated using the second of equation (1). One just 
needs the cross sections for 77 production of these pairs. Since they depend 
on the inverse of the square of the particle mass [63], the pair-production 
cross sections are much smaller in this case. The same applies to single 
meson production by 77 fusion. Low [97] showed that one can relate the 
particle production by two real photons (with energies lo\ and L02, respec- 
tively) to the particle's decay width, T 77 . Since both processes involve the 
same matrix elements, only the phase-space factors and polarization sum- 
mations are distinct. Low's formula is 

a{u!, uj 2 ) = 8tt 2 (2J + lfe 5(4^2 - M 2 ) , (28) 
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Figure 7: Feynman graphs for two- and three-photon fusion in ultra-peripheral 
collisions of relativistic heavy ions. 

where J, M, and Tm-*^ are the spin, mass and two-photon decay width 
of the meson, W is the cm. energy of the colliding photons [63]. The 
delta-function accounts for energy conservation. 

In ref. [5] the following equation was obtained for the production of 
mesons with mass M in HI colliders: 



where 5 = 0.681... Later [99] it was shown that a more detailed account of 
the space geometry of the two-photon collision is necessary, specially for the 
heavier mesons. R is a parameter which depends on the mass of the produced 
particle. If M is much smaller than the inverse of a typical nuclear radius, 
then R = 1/M, otherwise R is the nuclear radius. These choices reflect the 
uncertainty relation in the direction transverse to the beam, as explained in 
ref. [5]. Since spin 1 particles cannot couple to two real photons [98], one 
expects that only spin and spin 2 particles are produced. 

In the next section we show how one can calculate meson production 
using the concepts of QED extended to include bound quark states [96] . 




(29) 
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5.2 Two-photon fusion in heavy ion colliders 

In the laboratory frame the Fourier components of the classical electromag- 
netic field at a distance b/2 of nucleus 1 with charge Ze and velocity (3, is 
given by (in our notation q = (<?o, Qt, <?3), and q 3 = q z ) 



A^\q)=-M 2 Ze 5(q -f3q 3 



e iqt.b/2 

' <it+<&h 2 



and A 



(i) _ 



M 1} 



(30) 



For the field of nucleus 2, moving in the opposite direction, we replace (3 
by —(3 and b by — b in the equations above. Although (3 ~ 1 in relativistic 
colliders, it is important to keep them in the key places, as some of their 
combinations will lead to important 7 = (1 — [3 2 )~ 1 ^ 2 factors. 

The matrix element for the production of positronium is directly ob- 
tained from the corresponding matrix element for the production of a free 
pair (see fig. 7(a)), with the requirement that P + = P_ = P/2, where P is 
the momentum of the final bound state. That is 



M 



= M1 + M2 
= —ie 2 u(—) 



J (2vr) 4 * { 2 q) q 2 -MV4^ 1 2 + Q) 



+ 



4 (1) (y + s)- 
1 1 



(27T) 

' + M/2 



^ (2)( 2 



M 2 /4 



(31) 



where M is the positronium mass. 

The treatment of bound states in quantum field theory is a very complex 
subject (for reviews, see [100, 101]). In our case, we want to use the matrix 
element for free-pair production and relate the results for the production 
of a bound-pair. A common trick used in this situation is to convolute 
the matrix element given above with the bound-state wave function. One 
can show (see, e.g., [102]) that this is equivalent to the use of projection 
operators of the form 



u 



u 



*(0) 



2VM 
vP(O) 



2VM 



(P + M)i 7 5 
(P + M)i f\ 



and 



(32) 



where • • • is any matrix operator. The first equation applies to a spin 
(parapositronium) and the second to spin 1 (orthopositronium) parti- 
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cles, respectively. In these equations $(r) is the bound state wavefunc- 
tion calculated at the origin, and e* is the polarization vector, given by 
e* ±1 = (0, l/y/2, ±i/y/2, 0) and e* = (0, 0, 0, 1). 

Using eq. 32 in eq. 31, one gets for the parapositronium production 



M 



M 



P X I 



where 



/ 



d 2 qt qt 
<£ + Q 2 



(P t /2 + q t ) 2 + u, 2 / 7 2 (P t /2 - q t ) 2 + <4/i 



with 



Q 



M 2 

2 



j L j 2_ ~ 

^ 4 ^ 2 7 2 



M 2 if 

2 + 4 



(33) 



(34) 



(35) 



and 



E/P ~ Pz 



E/(3 + P z 



uj 2 = — — ^ and 4wiw 2 = M 2 +P 2 -P z /~f 2 ~ M 2 +P 2 

(36) 

where 22 = Po is the total positronium energy. 

We see that uj\ and a->2 play the role of the (real) photon energies. For 
real photons one expects Auiuj2 = E 2 , as in eq. 28. 

The two-photon fusion cross sections can be obtained by using 



da = J2 Jd 2 b\M(p)f 



d 3 P 



(2vr) 3 2E 



(37) 



The positronium wavefunction at the origin is very well known. It is given 
by ^(0) = M 3 a 3 /647r, where M is the positronium mass. 

Since the important impact parameters for the production of the positro- 
nium will be b > l/m e S> R, where R is the nuclear radius, the integral over 
impact parameter can start from 6=0. Thus, the integral over impact 
parameter in eq. 37 yields a delta function. Changing to the variables 



qi* 



-q* . 



qiz 



Pz/2 



eq. 37 becomes 



E 



da C(3) (Jo 



d?P 



7T 



M 2 



7 



Jb , 



Q2z 



Pz/2 + q z 



1 



(38) 
(39) 
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where 

M 



M [ 

Jb = / A 2 5 (qu + q 2 t - Pt) dqitdq 2 t , (40) 

7T J 



with 



qi t x qgf M 2 _ , 4Z 4 a 7 

g 2 g 2 M--^-^2 



A = ^x^ M2 _ , , , CO) = 1.202... , and „ _ ^ 



(41) 

We have included the zeta-function £(3) to take into account the production 
of the para-positronium in higher orbits, besides the production in the K- 
shell. 

The equation 39 is equal to the equation obtained in ref. [66] for the 
production of the parapositronium, using an alternative method. 

The integral in eq. 40 can be carried out analytically, reducing it to a 3- 
dimensional integral [96]. For RHIC, using 7 = 108 and Au + Au collisions, 
one finds a = 19 A mb. For the LHC , using 7 = 3000 and Pb+Pb collisions, 
one finds a = 116 mb. These are in good agreement with the results (Born 
cross sections) of ref. [66]. Notice however, that Coulomb corrections are 
very important, due to the low mass of the electron and the positron [66]. 
When Coulomb corrections are included to the Born cross sections, the final 
values decrease by as much as 43% for RHIC and 27% for LHC. For meson 
production one expects that these corrections ere less relevant. 

5.3 Production of C = even Mesons 

One can extend the calculation of the previous section to account for the 
production of mesons with spin J = and J = 2 by the two-photon fusion 
mechanism. The following procedure is to be adopted: 

1. Replace the electron-positron lines by quark- antiquarks in the diagram 
of figure 7(a). 

2. M in the following formulas will refer to the meson mass. 

3. Replace a 2 by a 2 (2 J + 1) 3J2iQfi where 3 accounts for the num- 
ber of colors, and Qi is the fractional quark charge. These two last 

2 

factors will cancel out when one expresses \£(0) in terms of T 77 , 
the decay-width of the meson. To understand how this is done, lets 
discuss the basics of the annihilation process of a positronium (see 
also ref. [104]). With probability a 2 the e~ can fluctuate and emit a 
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virtual photon with energy m e . The electron recoils and can travel 
up to a distance ~ l/m e (or time ~ m e ) to meet the positron and 
annihilate. This occurs when e~ and e + are both found close to- 
gether in a volume of size (l/m e ) 3 , i.e., with a probability given by 
| ^/ (0) | 2 / rri^ . Thus, the annihilation probability per unit time (decay 
width) is T ~ a 2 \^(0)\ 2 /m 2 . Angular momentum conservation and 
CP invariance does not allow the ortho-positronium to decay into an 
even number of photons [98] . The description of the annihilation pro- 
cess given above is thus only appropriate for the para-positronium. 
A detailed QED calculation yields an extra 4ir in the formula above. 
This yields T 77 ( 1 5o) = 8.03 x 10 9 s _1 , while the experimental value 
[107] is 7.99(11) x 10 9 in good agreement with the theory. For 
mesons, including the color and the charge factors, as described before, 
the relationship between 'I'(O) and T 77 arise due to the same reasons. 

One gets T 77 = 16vra 2 f(0) 2 /M 2 ■ 3£; Qf. 

2 

According to these arguments the connection between T 77 and ^(0) , 
extended to meson decays, should be valid for large quark masses so 
that l/m q <C V< r 2 >, where V < r 2 > is the mean size of the me- 
son. Thus, it should work well for, e.g. charmonium states, cc. In 
fact, Appelquist and Politzer [105] have generalized this derivation 
for the hadronic decay of heavy quark states, which besides other 
phase-space considerations amounts in changing a to a s , the strong 
coupling constant. This can be simply viewed as a way to get a con- 
straint on the wavefunction [^(O)) 2 (see, ref. [106]). One expects that 
these arguments are valid to zeroth order in quantum chromodynam- 
ics and in addition one should include relativistic corrections. But, as 
shown in [104], the inclusion of relativistic effects, summing diagrams 
to higher order in the perturbation series, is equivalent to solving the 
non-relativistic Schrodinger equation. 

4. Change the integration variable to qu and Q2t- 

5. Introduce form factors F(qu) and F(q2t) to account for the nuclear 
dimensions. This is a simple way to eliminate the integral over impact 
parameters and is justified 'a posteriori', i.e., when one compares the 
results with those from other methods. These form factors will impose 
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a cutoff in q± t and q2t , so that 

git , 92* ^ ^ < M , (42) 

where R is a typical nuclear size. Taking R = 6.5 fm, one gets 1/R ~ 
30 MeV. This is much smaller than the meson masses. As an outcome 
of this condition, one can replace Q 2 ~ M 2 /2 in eq. 35. 

According the procedures 1-5 one gets from eq. 39, 

da 16(2J + 1) Z A a 2 j_ 

~ V 2 F 77 1 

x /dcutdqa (qat x q 2 ,) 2 [ J, i(^)^(^)] 2 (43) 

Using eqs. 36 one has 

E = Ui + L0 2 , cji - u; 2 = P 2 , and ujiuj 2 = M 2 /4 (44) 

so that 

dP,= (l + ^W and E = ^ + M2 / 4 . (45) 

Thus, 

dcr (+) dA/2 7 (m) (+ ) 1 / \ / \ (AR , 

-. — = cr w — -f = a [ ^> — ni(wi)n 2 (a; 2 ) , (46) 

where 

<7« = 8vr 2 (2J+l)^l and n^) = - Z 2 a [ dg ^ . (47) 

We notice that n(u;) is the frequently used form of the equivalent photon 
number which enters eq. 1. Thus, eqs. 46 and 47 are the result one expects 
by using the equivalent photon method, i.e., by using eqs. 28 and 1. This is 
an important result, since it shows that the projection method to calculate 
the two-photon production of mesons works even for light quark masses 
(i.e., for 7r ). In this case there seems to be no justification for replacing the 
quark masses and momenta by half the meson masses and momenta, as we 
did for the derivation of eq. 43. This looks quite intriguing, but it is easy 
to see that the step 3 in our list of procedures adopted is solely dependent 



24 



on the meson mass, not on the quark masses, i.e., if they are constituent, 
sea quarks, etc. Moreover, the projection method eliminates the reference 
to quark masses in the momentum integrals. The condition 42 finishes the 
job, by eliminating the photon virtualities and yielding the same result one 
would get with the equivalent photon approximation. 

One can define a "two-photon equivalent number", A/2 7 (M 2 ), so that 



To calculate this integral one needs the equivalent photon numbers given 
by eq. 47. The simplest form factor one can use for this purpose is the 
'sharp-cutoff' model, which assumes that 

F{q 2 ) = 1 for q 2 < 1/R , and F{q 2 ) = , otherwise . (49) 

In this case, one gets for the differential cross section 



The spectrum possesses a characteristic 1/uj dependence, except for ui 3> 
j/R, when it decreases as 1/uj 5 . 

When the condition 7 3> MR is met, one can neglect the unity factors 
inside the logarithm in eq. 50, as well as the second term inside brackets. 
Then, doing the integration of 50 from u = M 2 R/4/j to uj = j/R, we get 
eq. 29. But eq. 50 is an improvement over eq. 29. Eq. 29 is often used 
in the literature, but it is only valid for 7 S> MR. This relation does not 
apply to, e.g., the Higgs boson production (M H o ~ 100 GeV), as considered 
in ref. [103, 108, 109]. 

For quantitative predictions one should use a more realistic form fac- 
tor. The Woods-Saxon distribution, with central density po, size R, and 
diffuseness a gives a good description of the densities of the nuclei involved 
in the calculation. For Au we use R = 6.38 fm, and a = 0.535 fm, with po 
normalized so that / d 3 rp(r) = 197. For Pb the appropriate numbers are 
6.63 fm, 0.549 fm, and 208, respectively [110]. In table 2 we show the cross 



a = 



aMjV 2y (M 2 ), where 
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Figure 8: Two-photon luminosities (see definition in eq. 51) at RHIC and 
LHC. Dashed lines include a geometric correction. 

sections for the production of C = even mesons at RHIC {Au + Au) and 
LHC (Pb + Pb) using the formalism described above. 

As pointed out in refs. [99, 111], one can improve the (classical) calcula- 
tion of the two-photon luminosities by introducing a geometrical factor (the 
B-function in ref. [99]), which affects the angular part of the integration 
over impact parameters. This factor takes care of the position where the 
meson is produced in the space surrounding the nuclei. In the approach 
described here the form factors also introduce a geometrical cutoff implying 
that the mesons cannot be produced inside the nuclei. However, it is not 
easy to compare both approaches directly, as one obtains a momentum rep- 
resentation of the amplitudes when one performs the integration over impact 
parameters to obtain eq. 43. But one can compare the effects of geometry 
in both cases by using equation 50. After performing the integral over to, 
one can rewrite it as 



where s = Auj\uj2 is the square of the center-of-mass energy of the two pho- 
tons, a(s) is given by eq. 28, and C(s) is the "photon-photon luminosity", 
given by 




(51) 
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Table 2: Cross sections for two-photon production of (C = even) mesons at 
RHIC (Au + Au) and at LHC (Pb + Pb). 



meson 


mass 


r 77 










a LHC 




[MeV] 


[kcV] 


[nb] 






[Mb] 


[mb] 


7T0 


134 


7.8 x 1CT 3 


99 


49 


2.8 


4940 


28 


V 


547 


0.46 j 


86 


12 


1.8 


1000 


16 


rf 


958 


4.2 


147 


5.1 


1.4 


746 


21 


/ 2 (1270) 


1275 


2.4 


179 


3.0 


1.2 


544 


22 


a 2 (1320) 


1318 


1.0 


67 


2.9 


1.1 


195 


8.2 


% 


2981 


7.5 


8.7 


0.38 


0.7 


3.3 


0.61 


XOc 


3415 


3.3 


2.6 


0.24 


0.63 


0.63 


0.16 


X2c 


3556 


0.8 


2.8 


0.21 


0.56 


0.59 


0.15 



/ 16 7 2 cj 2 \ 1 1 

x ln { 1 + ~1^W ) ~ TT^nJJW^l ■ (52) 

In figure 8 we compare the result obtained by eq. 52 and that of ref. [99] . 
The luminosities for RHIC (Au + Au) and for LHC (Pb + Pb) are presented. 
For RHIC the difference between the two results can reach 10% for very 
large meson masses (e.g. the Higgs), but one notices that for the LHC the 
two results are practically identical, the difference being of the order of 3%, 
or less, even for the Higgs production. Thus, the improved version of eq. 29, 
given by integrating eq. 50, is accurate enough to describe meson produc- 
tion by two-photon fusion. Other effects, like the interference between the 
electromagnetic and the strong interaction production mechanism in graz- 
ing collisions, must yield larger corrections to the (non-disruptive) meson 
production cross sections than a more elaborate description of geometrical 
effects. 

One might think that the calculation above can be extended to the three- 
photon fusion by using the equivalent photon approximation that, as we 
have seen in this section, works so well for C = even mesons. However, 
the introduction of a third photon leads to an additional integration, which 
implies that at least two of the exchanged photons cannot be treated as 
real ones. Nonetheless, the results of this section paves the way to the 
calculation of production of C = odd mesons. Although the use of the 
projection technique to systems composed of light quarks is questionable, 
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Figure 9: Feynman graphs for three-photon fusion in ultra-peripheral collisions of 
relativistic heavy ions. 

we have seen that it works, basically because of the relation 42, due to the 
inclusion of the nuclear form- factors. 



5.4 Three-photon Mechanism: Vector mesons 



Lets now consider the diagram of figure 7(b), appropriate for the fusion of 
three photons into a C = odd particle. According to the Feynman rules, the 
matrix element for it is given by 



Mr 



X ^ 2 \q 



d A q 



k) 



(2vr)4 
}6 + M/2 



(2tt)^ [ 2 



P 



k 2 



j+M/2 

q 2 -M 2 /4 



(53) 



There will be 12 diagrams like this. But the upper photon leg in diagram 
of fig. 7(b) can be treated as a real photon, meaning that the equivalent 
photon approximation is valid for this piece of the diagram. 

A straightforward calculation yields 



da 
dPz 



= 1024tt 



tf(0) (Za) 



1 



/ 



M 3 E 

2 



/ 



dqitqft [F(q!t)Y 



dq 2 t gg [Fjqlt)} 
[q 2 t + (2wx - u 2 ) 2 h 2 ] 



I 



(qft + ^h 2 ) 

dkt h F{k 2 ) 



1 2 



( fc 2 + (wi _ W2) 2 /7 2) 



(54) 



We now use the relationship between E and P z to oj\ and uj 2 and get rid 

2 

of the meson wavefunction at the origin. The wavefunction ^(0) cannot 
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be related to the 77 decay widths. But, vector mesons can decay into e + e~ 

pairs. These decay widths are very well known experimentally. Following a 

similar derivation as for the 77-decay the e + e~ decay-width of the vector 

2 

mesons can be shown [106] to be equal to T e+e - = Wira 2 ^(0) /3M 2 (3 • 

J2 i Q 2 ). Inserting these results in the above equation, the factor (3 • X)i Qf) 
will cancel out for the same reason as explained in section 3, and one gets 

H(M,lo) (55) 



da (_) n(uj) 



du u> 

where 

rr(-) — QP.TT- 

M 3 



* (_) = 96vr%F , (56) 



with n(to) given by 47 and 

H{M,u) = Z^M 2 I dq2t q2t [F{qlt)] 



1 



[4 + (M72^-„)7 7 2 ] 
dk t k t F(k 2 ) 



(k 2 + {M 2 /Auj -uj) 2 h 2 



(57) 



The above formulas should also be valid for the production of the ortho- 
positronium in ultra-peripheral collisions of relativistic heavy ions. For 
RHIC (Au + Au) one obtains a = 11.2 mb, while for the LHC (Pb + Pb) we 
get a = 35 mb. These numbers are also in good agreement with the results 
(in the Born approximation) given in ref. [66] . When one includes Coulomb 
corrections, as shown in ref. [66], the cross sections for orhto-positronium 
production is reduced by 40% for both RHIC and LHC. This is not consid- 
ered in the present approach as we are mainly interested in vector meson 
production for which this effect should be smaller. 

In table 3 we present the cross sections for the production of vector 
mesons by means of the three-photon fusion process are given. 

One sees that the cross sections for the production of vector mesons 
in ultra-peripheral collisions of relativistic heavy ions are small. They are 
not comparable to the production of vector mesons in central collisions. In 
principle, one would expect that the cross sections for three-photon produc- 
tion would scale as (Za) 3 , which is an extra Za factor compared to the 
two-photon fusion cross sections. However, the integral over the additional 
photon momentum decreases the cross section by several orders of magni- 
tude. 
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Table 3: Cross sections for three-photon production of vector (C = odd) 
mesons at RHIC (Au + Au) and at LHC (Pb + Pb). 



meson 


mass [MeV] 


r e+e - [kcv] 


a(-> [nb] 


a RHW [ nb ] 


a LHC [nb] 


p" 


770 


6.77 


1740 


137 


1801 




782 


0.60 


147 


13 


163 


J/V 


3097 


5.26 


21 


31 


423 


*' 


3686 


2.12 


5 


12 


155 



In the next section we discuss the possibility to access information on 
the gluon distribution in nuclei using the abundant virtual photons in PHIC 
[112]. 

6 One-photon Production Mechanism 

6.1 Vector Mesons and Gluon Distributions 

One of main predictions of QCD is the transition from the confined/chirally 
broken phase to the deconfined/chirally symmetric state of quasi-free quarks 
and gluons, the so-called quark-gluon plasma (QGP). Recently the heavy ion 
collisions have provided strong evidence for the formation of a QGP [113], 
with the first results of RHIC marking the beginning of a collider era in 
the experiments with relativistic heavy ions, as well as the era of detailed 
studies of the characteristics of the QGP. Currently, distinct models asso- 
ciated to different assumptions describe reasonably the experimental data 
[114], with the main uncertainty present in these analysis directly connected 
with the poor knowledge of the initial conditions of the heavy ion collisions. 
Theoretically, the early evolution of these nuclear collisions is governed by 
the dominant role of gluons [115], due to their large interaction probabil- 
ity and the large gluonic component in the initial nuclear wave functions. 
This leads to a "hot gluon scenario" , in which the large number of initially 
produced energetic partons create a high temperature, high density plasma 
of predominantly hot gluons and a considerably number of quarks. Such 
extreme condition is expected to significantly influence QGP signals and 
should modify the hard probes produced at early times of the heavy ion 
collision. Consequently, a systematic measurement of the nuclear gluon dis- 
tribution is of fundamental interest in understanding of the parton structure 
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of nuclei and to determine the initial conditions of the QGP. Other impor- 
tant motivation for the determination of the nuclear gluon distributions is 
that the high density effects expected to occur in the high energy limit of 
QCD should be manifest in the modification of the gluon dynamics. 

At small x and/or large A we expect the transition of the regime de- 
scribed by the linear dynamics (DGLAP, BFKL) (for a review, see e.g. ref. 
[116]), where only the parton emissions are considered, to a new regime 
where the physical process of recombination of partons becomes important 
in the parton cascade and the evolution is given by a nonlinear evolution 
equation. In this regime a Color Glass Condensate (CGC) is expected to be 
formed [117], being characterized by the limitation on the maximum phase- 
space parton density that can be reached in the hadron/nuclear wavefunc- 
tion (parton saturation) and very high values of the QCD field strength 
Ffiv ~ l/\/S7 [118]. In this case, the number of gluons per unit phase space 
volume practically saturates and at large densities grows only very slowly 
(logarithmically) as a function of the energy [119]. This implies a large mod- 
ification of the gluon distribution if compared with the predictions of the 
linear dynamics, which is amplified in nuclear processes [120, 121, 122, 123]. 

Other medium effects are also expected to be present in the nuclear 
gluon distribution at large values of x: the antishadowing (0.1 < x < 0.3), 
the EMC effect (0.3 < x < 0.7) and the Fermi motion (x > 0.7) [124, 125]. 
The presence of these effects is induced from the experimental data for the 
nuclear structure function which determines the behavior of the nuclear 
quark distributions and the use of the momentum sum rule as constraint. 
Experimentally, the behavior of the nuclear gluon distribution is indirectly 
determined in the lepton-nucleus processes in a small kinematic range of 
the fixed target experiments, with the behavior at small x (high energy) 
completely undefined. This situation should be improved in the future with 
the electron-nucleus colliders at HERA and RHIC [126, 127], which proba- 
bly could determine whether parton distributions saturate. However, until 
these colliders become reality we need to consider alternative searches in the 
current accelerators which allow us to constraint the nuclear gluon distribu- 
tion. In this section we analyze the possibility of using peripheral heavy ion 
collisions as a photonuclear collider and therefore to determine the behavior 
of the gluon distribution. 

A photon stemming from the electromagnetic field of one of the two 
colliding nuclei can penetrate into the other nucleus and interact with one 
or more of its hadrons, giving rise to photon-nucleus collisions to an energy 
region hitherto unexplored experimentally. For example, the interaction of 
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Figure 10: Energy dependence of the photoproduction of heavy quarks for 
distinct nuclear gluon distributions (A = 208). 

quasireal photons with protons has been studied extensively at the electron- 
proton collider at HERA, with yfs = 300 GeV. The obtained 7p center of 
mass energies extends up to W 7P « 200 GeV, an order of magnitude larger 
than those reached by fixed target experiments. Due to the larger number 
of photons coming from one of the colliding nuclei in heavy ion collisions 
similar and more detailed studies will be possible in these collisions, with 
W*yN reaching 950 GeV for the Large Hadron Collider (LHC) operating in 
its heavy ion mode. 

When a very hard photon from one equivalent swarm of photons pen- 
etrates the other nucleus it is able to resolve the partonic structure of the 
nucleus and to interact with the quarks and gluons. One of the basic process 
which can occur in the high energy limit is the photon-gluon fusion leading 
to the production of a quark pair. The main characteristic of this process is 
that the cross section is directly proportional to the nuclear gluon distribu- 
tion. The analysis of this process in peripheral heavy ion collisions has been 
proposed many years ago [128] and improved in the Refs. [129, 130] (for a 
review, see ref. [131]). Here we reanalyze the charm photoproduction as a 
way to estimate the medium effects in xGa in the full kinematic region. We 
consider as input three distinct parameterizations of the nuclear gluon distri- 
bution. First, we consider the possibility that the nuclear gluon distribution 
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is not modified by medium effects, i.e., xGa(x, Q 2 ) = Ax xGn(x, Q 2 ), with 
the nucleon gluon distribution (xGn) given by the GRV parameterization 
[132]. Moreover, we consider that xGa{x,Q 2 ) = Rg x xGn{x,Q 2 ), where 
Rg parameterizes the medium effects as proposed by Eskola, Kolhinen and 
Salgado (EKS) [133]. The main shortcoming of these parameterizations is 
that these are based on the DGLAP evolution equations which are not valid 
in the small x regime, where the parton saturation effects should be consid- 
ered. In order to analyze the sensitivity of peripheral heavy ion collisions to 
these effects we consider as input the parameterization proposed by Ayala 
Filho and Goncalves (AG) [134] which improves the EKS parameterization 
to include the high density effects. 

One shortcoming of the analysis of photoproduction of heavy quarks in 
peripheral heavy ion collisions to constraint the nuclear gluon distribution 
is the linear dependence of the cross section with this distribution. This 
implies that only experimental data with a large statistics and small error 
will allow to discriminate the medium effects in the nuclear gluon distribu- 
tion. Consequently, it is very important to analyze other possible processes 
which have a stronger dependence in xGa- Here we propose the study of the 
elastic photoproduction of vector mesons in peripheral heavy ion collisions 
as a probe of the behavior of the nuclear gluon distribution. This process 
has been largely studied in ep collisions at HERA, with the perturbative 
QCD predictions describing successfully the experimental data [116], con- 
sidering a quadratic dependence of the cross section with the nucleon gluon 
distribution. 

6.2 Photoproduction of Heavy Quarks 

At high energies the dominant process occurring when the photon probes the 
structure of the nucleus is the photon-gluon fusion producing a quark pair. 
For heavy quarks the photoproduction can be described using perturbative 
QCD, with the cross section given in terms of the convolution between the 
elementary cross section for the sub-process — > QQ and the probability 
of finding a gluon inside the nucleus, i.e., the nuclear gluon distribution. 
Basically, the cross section for cc photoproduction is given by 



where da C c/dM C c is calculable perturbatively, is the invariant mass of 
the cc pair with M 2 Z = s = xs, s is the squared CM energy of the jA 




(58) 



33 




2 6 10 14 18 22 26 30 
M (GeV) 



Figure 11: Differential cross section for cc production versus the invariant 
mass M = M^. The photon- gluon fusion in the heavy- ion collision system 
208 Pb + 208 Pb at LHC energy is considered. 



system, g^(x,fj,) is the gluon density inside the nuclear medium, fj, is the 



factorization scale (fi = J M^) , and m c is the charm quark mass (we assume 
that m c = 1.45 GeV). Moreover, the differential cross section is [135] 



dMrr ~ M, 



cc 



(1 + 6 + le 2 )ln( | + ^ZI )- (i + e )vT^] , (59) 



where e c is the charm charge and e = 4m 2 /M 2 -. From the above expression, 
we verify that the cross section is directly proportional to the nuclear gluon 
distribution, which implies the possibility to constraint its behavior from 
experimental results for photoproduction of heavy quarks. 

Throughout this section we use the Born expression for the elementary 
photon-gluon cross section [eq. (59)]. QCD corrections to the Born cross 
section will not be considered here, although these corrections modify the 
normalization of the cross section by a factor of two. This is justified by the 
fact that we are interested in the relative difference between the predictions 
of the distinct nuclear gluon distributions, which should be not modified by 
the next-to-leading-order corrections. 

In Fig. 10 we present the energy dependence of the photoproduction 
cross section. We focus our discussion on charm photoproduction instead 
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of bottom photoproduction, since in this process smaller values of x are 
probed. One verifies that different nuclear gluon distributions imply distinct 
behaviors for the cross section, with the difference between the predictions 
increasing with the energy. This is associated to the fact that at high ener- 
gies we are probing the small x behavior of xGa, since x oc M C c/s, where 
M C c is the invariant mass of the photon-gluon system. Currently, only the 
region of small center of mass energy has been analyzed by the fixed tar- 
get electron- nucleus experiments, not allowing a good constraint on medium 
effects present in the nuclear gluon distribution. Such situation should be 
improved in the future with electron-nucleus colliders at HERA and RHIC 
[124, 127]. 

Another possibility to study photoproduction of heavy quarks at large 
center of mass energies is to consider peripheral heavy ion collisions [128, 
129, 130] . In this process the large number of photons coming from one of the 
colliding nuclei in heavy ion collisions will allow to study photoproduction, 
with W 7 jv reaching 950 GeV for the LHC. To determine the photoproduction 
of heavy quarks in peripheral heavy ion collisions the elementary photon- 
gluon cross section has to be convoluted with the photon energy distribution 
and the gluon distribution inside the nucleus: 

a {A A -► XXQQ) = n(u) ® ® xG A (x, Q 2 ) . (60) 

It is interesting to determine the values of x which will be probe in periph- 
eral heavy ion collisions. The Bjorken x variable is given by x = (M/2p)e~ y , 
where M is the invariant mass of the photon-gluon system and y the center 
of momentum rapidity. For Pb + Pb collisions at LHC energies the nucleon 
momentum is equal to p = 3000 GeV; hence x = (M/6000 GeV)e~ y . There- 
fore, the region of small mass and large rapidities probes directly the small x 
behavior of the nuclear gluon distribution. This demonstrates that periph- 
eral heavy ion collisions at LHC represents a very good tool to determine the 
behavior of the gluon distribution in a nuclear medium, and in particular 
the low x regime. Conversely, the region of large mass and small rapidi- 
ties is directly associated to the region where the EMC and antishadowing 
effects are expected to be present. Similarly, for RHIC energies (p = 100 
GeV) the cross section will probe the region of medium and large values of 
x (x > 10~ 2 ). For this kinematical region, the EKS and AG parameteriza- 
tions are identical, which implies that the photoproduction of heavy quarks 
in peripheral heavy ion collisions at RHIC does not allows to constraint the 
high density effects. However, the study of this process at RHIC will be 
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Figure 12: Rapidity distribution for the photoproduction of charm quarks 
in 208 Pb + 208 Pb collisions at LHC. 

very interesting to determine the presence or not of the antishadowing and 
EMC effect in the nuclear gluon distribution. 

Fig. 11 shows the mass distribution in photoproduction of charm quarks 
in peripheral heavy ion collisions. We can see that the main difference 
between the predictions occur at small values of M, which is associated 
to the small x region. Basically, we have that the predictions of the EKS 
parameterization are a factor 1.25 larger than the AG prediction in this 
region, while the prediction of Gliick, Reya and Vogt (GRV) is a factor 
2.4 larger. This result is consistent with the fact that the main differences 
between the parameterizations of the nuclear gluon distribution occur at 
small x [134]. The difference between the predictions diminishes with the 
growth of the invariant mass, which implies that this distribution is not a 
good quantity to estimate the nuclear effects for medium and large x. 

A better distribution to discriminate the behavior of the nuclear gluon 
distribution is the rapidity distribution, which is directly associated to the 
Bjorken x variable, as discussed above. The rapidity distribution is calcu- 
lated considering that da/dy = ujda/duj. In Fig. 12 we show the rapidity 
distribution considering the three parameterizations of xGa as input. We 
have that the region of small rapidities probes the region of large x, while for 
large y we directly discriminate the different predictions of xGa for small x. 
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Figure 13: Rapidity distribution for the photoproduction of charm quarks 
in 208 Pb + 208 Pb collisions at LHC with (GRV+RES) and without (GRV) 
the inclusion of the resolved contribution. 



These results are coherent with this picture: as at large x the EKS and AG 
predictions are identical, this region allows to estimate Rq = xGa/(AxGn) 
in the region of the antishadowing and EMC effects; at large y the large 
difference between the parameterizations implies large modifications in the 
rapidity distribution, which should allow a clean experimental analysis. 

A comment is in order here. In hard photon-hadron interactions the 
photon can behave as a pointlike particle in the so-called direct photon pro- 
cesses or it can act as a source of partons, which then scatter against partons 
in the hadron, in the resolved photon processes (for a recent review see ref. 
[136]). Resolved interactions stem from the photon fluctuation to a quark- 
antiquark state or a more complex partonic state, which are embedded in 
the definition of the photon structure functions. Recently, the process of 
jet production in photoproduction has been search of studies of the par- 
tonic structure of the photon (See e.g. [116]), and the contribution of the 
resolved photon for the photoproduction of charm has been estimated [137]. 
Basically, these studies shown that the partonic structure of the photon is 
particularly important in some kinematic regions (for example, the region 
of large transverse momentum of the charm pair [137]). Consequently, it is 
important to analyze if this contribution modifies, for instance, the results 
for the rapidity distribution, which is strongly dependent of nuclear gluon 
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distribution. To leading order, beyond the process of photon-gluon fusion 
considered above, charm production can occur also in resolved photon inter- 
actions, mainly through the process gg — > cc. Therefore, one needs to add 
in eq. (60) the resolved contribution given by 

a res (AA — ► XXQQ) = n(u) <8> x 1 G 1 {x 1 , Q 2 ) ® v gg ^QQ ®xG A (x,Q 2 ) ,(61) 

where x 7 denote the fraction of the photon momentum carried by its gluon 
component x 7 G 7 and a g g ^QQ the heavy quark production cross section first 
calculated in ref. [138]. Since the resolved contribution should be the same 
for the three nuclear parton distributions, Fig. 13 only shows the results 
obtained using the GRV parameterization for the nucleon. For the photon 
distribution we use the GRV photon parameterization [139], which predicts 
a strong growth of the photon gluon distribution at small x 1 . We can see 
that though this contribution is important in the photoproduction of heavy 
quarks, as shown in ref. [137], it is small in the rapidity distribution of 
charm quarks produced in peripheral heavy ion collisions. Therefore, the 
inclusion of the resolved component of the photon does not is not relevant 
for the use of this process as a probe of the nuclear gluon distribution. 

It is important to salient that the potentiality of the photoproduction of 
quarks to probe the high density effects have been recently emphasized in ref. 
[140], where the color glass condensate formalism was used to estimate the 
cross section and transverse momentum spectrum. The authors have verified 
that the cross section is sensitive to the saturation scale which characterizes 
the colored glass. The results abovecor rob orate the conclusion that this 
process is sensitive to the high density effects and the agreement between 
the predictions is expected in the kinematic region in which the transverse 
momentum of the pair kt is larger than the saturation scale Q s . For k t < Q s 
the collinear factorization to calculate the cross sections must be generalized, 
similarly to ref. [140]. 

6.3 Elastic Photoproduction of Vector Mesons 

The production of vector mesons at HERA has become a rich field of ex- 
perimental and theoretical research [for a recent review see, e.g. ref. [141]], 
mainly related with the question of whether perturbative QCD (pQCD) can 
provide an accurate description of the elastic photoproduction processes. 
At high energies the elastic photoproduction of vector mesons is a two-stage 
process: at first the photon fluctuates into the vector meson which then 
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Figure 14: Energy dependence of the elastic photoproduction of J/\& for 
distinct nuclear gluon distributions (A = 208). 

interacts with the target. For light vector mesons the latter process oc- 
curs similarly to the soft hadron-hadron interactions and can be interpreted 
within Regge phenomenology [142]. However, at large mass of the vector 
meson, for instance the mass of the J/^> meson, the process is hard and 
pQCD can be applied [143]. In this case the lifetime of the quark-antiquark 
fluctuation is large compared with the typical interaction time scale and the 
formation of the vector meson only occur after the interaction with the tar- 
get. In pQCD the interaction of the qq pair is described by the exchange of 
a color singlet system of gluons (two gluons to leading order) and, contrary 
to the Regge approach, a steep rise of the vector meson cross section is pre- 
dicted, driven by the gluon distribution in the proton. Measurements of the 
elastic photoproduction of J/^ mesons in ep processes has been obtained 
by the HI and ZEUS collaborations for values of center of mass energy be- 
low 300 GeV, demonstrating the steep rise of the cross section predicted by 
pQCD. This result motivates the extension of the pQCD approach used in 
electron-proton collisions to photonuclear processes. 

The procedure for calculating the forward differential cross section for 
photoproduction of a heavy vector meson in the color dipole approximation 
is straightforward. The calculation was performed some years ago to leading 
logarithmic approximation, assuming the produced vector-meson quarko- 
nium system to be nonrelativistic [143] and improved in distinct aspects 
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Figure 15: The rapidity distribution for elastic photoproduction of J/^f at 
LHC considering distinct nuclear gluon distributions (A = 208). 

[144]. To lowest order the 7 A — * J/^fA amplitude can be factorized into 
the product of the 7 — > cc transition, the scattering of the cc system on the 
nucleus via (colorless) two-gluon exchange, and finally the formation of the 
J/ty from the outgoing cc pair. The heavy meson mass Mjw ensures that 
pQCD can be applied to photoproduction. The contribution of pQCD to 
the imaginary part of the t = differential cross section of photoproduction 
of heavy vector mesons is given by [143] 

J t 1*=° = ^^ =^x[xG A (x,Q)] , (62) 

where xGa is the nuclear gluon distribution, x = AQ 2 /W 2 with W the center 
of mass energy and Q 2 = M 2 ^/4. Moreover, T ee is the leptonic decay width 
of the vector meson. The total cross section is obtained by integrating over 
the momentum transfer t, 

c{lA _ Jm) = MiA-.j/*A) ^ r- df |F(()|2 m 

J train 

where t m i n = (M 2 j^/2uj) 2 and F(t) = J d 3 r p(r) exp( iq • r) is the nuclear 
form factor for the matter distribution. 
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Table 4: The total cross section a(AA — > AAJ/^f) for different nuclear gluon 
distributions. Results for LHC. 



Gluon Distribution 


LHC 


GRV 


6.584 mb 


EKS 


2.452 mb 


AG 


0.893 mb 



A comment is in order here. Although some improvements of the ex- 
pression (62) have been proposed in the literature [144], these modifications 
do not change the quadratic dependence on xGa- 

The main characteristic of the elastic photoproduction of vector mesons 
is the quadratic dependence on the gluon distribution, which makes it an 
excellent probe of the behavior of this distribution. In Fig. 14 we show the 
energy dependence of the differential cross section [eq. (62)], considering 
the three distinct nuclear gluon distributions discussed above. One obtains 
larger differences between the predictions than obtained in photoproduction 
of heavy quarks, mainly at large values of energy. This result motivates 
experimental analysis of elastic J/\I/ photoproduction in photonuclear pro- 
cesses at high energies. Although the future electron ion colliders (HERA-A 
and eRHIC) should probe this kinematic region, here we show that this can 
also be done with peripheral heavy ion collisions. 

Following similar steps used in photoproduction of heavy quarks, photons 
coming from one of the colliding nuclei may interact with the other. For 
elastic photoproduction of J/\l/ one can consider that this photon decays into 
a cc pair which interacts with the nuclei by the two gluon exchange. After the 
interaction, this pair becomes the heavy quarkonium state. Consequently, 
the total cross section for J/\I/ production in peripheral heavy ion collisions 
is obtained by integrating the photonuclear cross section [eq. (63)] over the 
photon spectrum, resulting 

«(AA AAJM = \*± n( ^M-> ^) |fao r rf(|F(t)|2 (64) 
J u at J tmin 

Table 4 presents the total cross section considering as input the distinct 
nuclear gluon distributions and LHC energies. Although these numbers 
will be modified by the inclusion of higher order corrections for the cross 
section [144], the difference between the predictions should not be altered. 
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Figure 16: Rapidity behavior of the ratio between the EKS and AG predic- 
tions for photoproduction of charm quarks and elastic photoproduction of 



One verifies that due to the large number of equivalent photons and the 
large center of mass energies of the photon-nucleus system, the cross section 
for this process is large, which allows an experimental verification of these 
predictions. Also, in peripheral heavy ion collisions the multiplicity is small 
what might simplify the experimental analysis. Moreover, the difference 
between the predictions is significant. For RHIC energies, the cross section 
for this process is small and probably an experimental determination will 
be very hard. 

In fig. 15 the rapidity distribution for J/\& production at LHC energies 
is shown. In this case, the final state rapidity is determined by 
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(65) 



Similarly to the heavy quark photoproduction, the large y region probes the 
behavior of xGa at small x, while the region of small y probes medium val- 
ues of x. One concludes that the rapidity distribution for elastic production 
of J/ty at RHIC allows to discriminate between the GRV and EKS predic- 
tion, with the AG prediction being almost identical to the latter. For LHC 
one finds a large difference between the distributions, mainly in magnitude, 
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which will allow to estimate the magnitude of the EMC, antishadowing and 
high density effects. 

Finally, fig. 16 compares the photoproduction of heavy quarks and the 
elastic photoproduction of J/^f in peripheral heavy ion collisions as a possi- 
ble process to constraint the behavior of the nuclear gluon distribution. The 
rapidity distribution of the ratio 

RIEKS/AG] = ^[EKS]/^[AG] , (66) 
dy dy 

is shown where we consider the EKS and AG parameterizations as inputs 
of the rapidity distributions. One confirms that the analysis of the elastic 
photoproduction of J/\I/ at medium and large rapidities is a potential process 
to determinate the presence and estimate the magnitude of the high density 
effects. 



7 Final Remarks 

The possibility to produce a Higgs boson via 77 fusion was suggested in ref. 
[103, 108]. The cross sections for LHC are of order of 1 nanobarn, about 
the same as for gluon-gluon fusion. But, the two-photon processes can also 
produce bb pairs which create a large background for detecting the Higgs 
boson. A good review of these topics was presented in ref. [145]. 

The excitation of a hadron in the field of a nucleus is another useful tool 
to study the properties of hadrons. It has been used for example to obtain 
the lifetime of the £0 particle by measuring the (Ml) excitation cross section 
for the process 7 + A — > So [146]. The vertex 7 — > 3ir has been investigated 
[147] in the reaction of pion pair production by pions in the nuclear Coulomb 
field: 7r~ + Z — > 7r~ + 7T° + Z. Also, the 7r~ polarizability has been studied 
in the reaction 7r~+Z^ir~+"/ + Z [148]. Other unexplored possibilities 
includes the excitation of a nucleon to a A-particle in the field of a heavy 
nucleus in order to disentangle the Ml and E2 parts of the excitation. 

As for meson-production in PHIC there are several planned experiments 
at RHIC, as well as for the future LHC [149]. These machines were de- 
signed to study hadronic processes. But, as have shown in this brief review, 
they can also be used for studying very interesting phenomena induced in 
peripheral collisions. 
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